D-A032  750 


unclassified 

I o f I 

A0A032750 


SOUTHEASTERN  MASSACHUSETTS  UNIV  NORTH  DARTMOUTH  DEPT  — ETC  F/G  20/6 
A MODIFIED  GRADIENT  TECHNIQUE  FOR  IMAGE  ANALYSIS. (U) 

SEP  76  C H CHEN  AF-AFOSR-2951-76 


EE-76-5 


AFOSR-TR-76-1205 


TECHNICAL  REPORT  SERIES  IN  INFORMATION  SCIENCES 
(Dr.  C.  H.  Chen,  Principal  Investigator) 


SOUTHEASTERN  MASSACHUSETTS 

UNIVERSITY 

ELECTRICAL  ENGINEERING  DEPARTMENT 


NORTH  DARTMOUTH,  MASS.  02747  U.S.A. 


I 


AIK  roses  0F7TC2  0 T SCT2STIFIC  RESEARCH  (ATSC) 
NOTICE  07  TRA35KITTA1.  TO  DOC  \ 

T'nia  technical  report  has  been  reviewed  and  is 
approved  for  public  release  IAW  AIR  190-12  (7b)  • 
Distribution  is  unlioited. 

A*  D*  BL05E 

Technical  Infornatlea  Officer 


v 


OT I 


Grant  AFOSR  76-2953 
TP  No.  EE-76-5 

September  15,  1976 


A MODIFIED  GF.ADIENT 
TECHNIQUE  FOR  IMAGE  ANALYSIS 


t>y 

C.  H.  Chen 

Electrical  Engineering  Department 
Southeastern  Massachusetts  University 
North  Dartmouth,  Massachusetts  02fhf 


Abstract 


The  modified  gradient  operation  considered  in  this  report  takes  the 
product  of  four  conventional  gradient  operations  in  four  different  directions. 
Such  modified  procedure  is  shown  by  extensive  computer  results  to  provide 
great  improvement  in  noise  reduction  and  boundary  extraction  over  the 
conventional  gradient  operation.  Furthermore,  the  procedure  is  computationally 
simple  and  insensitive  to  rotation  of  images.  Computer  results  of  modified 
gradient  operation  followed  by  image  compression  are  also  presented. 
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A Modified  Gradient  Technique  for  Image  Analysis 

C.  fi.  Chen 

The  gradient  operation  may  be  considered  as  an  approximation  to  the  spatial 
differentiation  of  grayscale  pictures.  The  resulting  gradient  picture,  however, 
has  a number  of  isolated  or  bursty  noise  elements,  as  well  as  disconnected  object 
boundaries.  The  modified  gradient  proposed  in  this  report  places  a greater 
emphasis  on  boundary  points,  which  usually  have  larger  gradient  values,  and  greatly 
suppresses  the  noises  resulting  from  differentiation.  To  comnress  the  image 
spatially,  it  is  more  suitable  to  use  the  modified  binary  gradient  picture 
because  more  boundary  information  is  available.  Compression  provides  further 
noise  suppression  in  addition  to  size  reduction. 

A variety  of  digital  gradient  operators  can  be  used  for  spatial  differentiation. 
For  the  -point  array 
A B 
C D 

the  gradient  may  be  defined  as  one  of  the  following  expressions, 

|<|A  - D|  + |B  - C|)  (la) 

|(|A  - B|  + | A - C| ) (lb) 

|(/(A  - B)2  + (A  - C)2  (lc) 

max( | C - D|,  |C  - A|)  (id) 

max( Ja  - D| , |B  - C| ) (le) 

Equation  (la)  is  the  most  commonly  used  conventional  gradient  operation. 

For  the  9-point  array  given  by: 

ABC 
D E F 


G H I 


-2- 

some  gradient  operators  are  defined  as 

max( |E  - A| , |E  - B| , |E  - c| , |E  - D| , | E - F | , 

|E  - G|,  |E  - H|,  |E  - l| ) (2a) 

max{ |(A  + B + C)  - (G  + H+  l)|,  i(A  + D + G)  - (C  + F + I)|}/3  (2b) 

/(D  - F)2  + (B  - H)2  + I (A  - I)2  + ( C - G)2  (2c) 

The  actual  performances  of  the  operations  given  by  Eqs.  (l)  and  (2)  are  very 
similar.  Theoretically  speaking,  an  infinite  array  would  be  needed  to  acliieve 
perfect  differentiation.  This  is  obviously  inroossible.  Consider  now  a l6-po.int 
array, 

A B C D 
E F G H 
I J K L 
M II  0 P 

The  proposed  modified  gradient  is  defined  as 


where 

a - |F  - K|  + |J  - g| 

b = J A — P | + |M  - D| 

c = | B - 0|  + 1 1 - II | 

a ® |c  - n|  + |e  - l| 

Each  point  of  the  array  is  used  once  and  only  once.  Equation  (3)  consists 
of  the  product  of  four  conventional  gradient  operation  (Eq.  la)  in  four  different 
directions.  (Unless  the  arrays  are  uniformly  speced  horizontally  and  vertically, 
a and  b do  not  coincide  in  directions).  Thus  the  modified  gradient  is  not 
sensitive  to  the  rotation  of  images.  Points  along  the  object  boundary  of  the 
original  picture  will  be  greatly  emphasized  because  the  product  of  large  gradient 
values  is  large.  Isolated  noise  points  will  be  reduced  because  at  least  one  of 
the  four  quantities  a,  b,  c,  d may  become  vary  small  and  the  product  would  be 


i. 


small  also. 


V 


BP 
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In  the  computer  study,  the-  display  in  provided  by  XY -plotter  which  hr.3  only 
two  levels:  pen-down  for  a black  point,  and  pen  up  for  blank.  A threshold  is 
set  such  taht  the  pen  is  down  when  the  threshold  is  exceeded.  Thus  only  binary 
picture  is  available.  In  Figs.  1-17,  the  gradient  picture  refers  to  a binary 
picture  computed  from  the  conventional  gradient  operation  given  by  Eq.  la;  and 
the  modified  gradient  picture  refers  to  a binary  picture  computed  from  Eq.  (3). 
Figures  1-9  show  a comparison  between  the  gradient  pictures  and  the  modified 
gradient  pictures.  For  a detailed  discussion  of  the  images  studied,  the  readers 
should  refer  to  the  report  by  this  author  entitled,  "Theory  and  Applications  of 
Image  Pattern  Recognition",  TR  EE-75-3,  April  1975,  prepared  for  Grant  AFOCR  71-2119- 

The  modified  gradient  obviously  provides  significant  improvement  over  the 
conventional  gradient  in  terms  of  noise  suppression  and  object  extraction.  Hie 
thickened  object  boundaries  may  be  thinned  by  thinning  algorithms.  Hie  smaller 
number  of  gaps  remaining  nay  be  filled  in  by  boundary  tracking  or  following 
algorithms.  The  required  amount  of  computation  which  is  slightly  more  than  four 
times  of  the  conventional  gradient  is  quite  small  as  compared  with  most  image 
analysis  algorithms. 

As  a further  step  in  image  analysis,  spatial  compression  from  the  modified 
binary  gradient  picture  may  be  obtained  by  reducing  every  9 points  (a  3 by  3 array) 
into  one  point.  A threshold  of  3 is  chosen  such  that  when  the  number  of  l's 
exceeds  or  equals  to  3 in  the  9 points,  the  9 points  are  replaced  by  a single 
point  of  1 (black):  otherwise  the  new  point  is  0 (blank).  Figures  10-17  show 
the  full  pictures  of  images  obtained  by  modified  gradient  operation  followed  by 
9 to  1 image  compression.  The  complete  operation  very  much  preserves  the  original 
boundary  information  provided  by  modified  gradient  while  the  additional  noise 
suppression  is  evident.  It  i3  noted  that  the  XY  plotter  plots  one  line  in  every 
12  lines.  Thus  there  are  actually  3^0  horizontal  lines  each  in  Figs.  10-17.  If 


A 


all  line3  are  displayed  the  objects  in  pictures  would  appear  to  be  more  continuous. 


From  the  compressed  pi ctures  we  notice  that  .'anther  improvement  should  he 
made  so  that  the  object  in  each  picture  cm  be  more  clearly  distinguished  from 
the  background.  One  other  problem  is  the  threshold  selection  to  generate  binary 
gradient  or  binary  modified  gradient  picture.  Figure  10  is  a typical  plot  of 


the  histograms  of  gradient  and  modified  gradient  pictures.  The  modified  gradient 
histogram  has  a sharper  peak  but' longer  tail.  The  tail  portion  corresponds  to  the 
object  boundary.  Thus  for  the  same  threshold  value,  the  area  above  threshold  is 
larger  in  the  modified  gradient  histogram  than  in  the  gradient  histogram.  If*  we 
lower  the  threshold  in  gradient  picture  to  get  the  same  number  of  points,  there 
would  be  considerable  increase  of  noise.  'This  explains  why  the  modified  gradient 
performs  better.  Hie  threshold  may  be  determined  by  specifying  a percentile 
(say*  15?)  of  points  above  threshold  or  simply  by  trial  and  error  on  the  scope  of 
the  keyboard  or  the  XY  plotter  itself. 
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The  modified  gradient  picture  corresponding 
to  Fig.  la  also  with  threshold  30. 


Fig.  3b  The  modified  gradient  picture  correspondirr; 
to  Fig  . 3a  also  witn  threshold  3 'j  . 


Fig.  4a  Tne  gradient  picture  of  the  target  area 
in  rile  7 with  threshold  30. 
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rig.  4b  The  modified  gradient!,  picture  corresponding 
to  Fig  ,4a  also  with  threshold  30  . 
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Fig.  7fc  The  modified  gradient  picture  corresponding 
to  Fig.  7a  also  with  threshold  50, 
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I- iq . 9 (a)  The  giadient  picture  and  ( 1 » ) the  modified  gradient  picture  of  the 

target  area  in  File  1()  with  threshold  J‘>  in  each  case. 


12  Full  picture  of  File  5:  modified  gradient  with  threshold  30 
followed  by  9 to  1 compression  of  threshold  3, 


13  Full  picture  of  File  7:  modified  gradient  with  threshold  30 
foil  »-ed  Ly  9 to  1 compression  of  threshold  3. 
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Full  picture  of  Fue  9:  nodified  gradient  with  threshold  45 
f . .lowed  by  9 to  1 compression  of  threohold  3. 
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16  Full  picture  of  File  11:  modified  gradient  with  threshold  30 
followed  by  9 to  1 compression  of  threshold  3. 
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Fig.  17  Full  picture  of  File  13:  modified  gradient  with  tnreshold  50 

followed  by  9 to  1 compression  of  threshold  3. 
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